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In classical general relativity, the generic approach to the initial singularity is very complicated 
as exemplified by the chaos of the Bianchi IX model which displays the generic local evolution close 
to a singularity. Quantum gravity effects can potentially change the behavior and lead to a simpler 
initial state. This is verified here in the context of loop quantum gravity, using methods of loop 
quantum cosmology: the chaotic behavior stops once quantum effects become important. This is 
consistent with the discrete structure of space predicted by loop quantum gravity. 
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According to the celebrated singularity theorems of 
classical general relativity, the backward evolution of an 
expanding universe leads to a singular state where the 
classical theory ceases to apply. An extensive analysis of 
the approach to the singularity, in the general context of 
inhomogeneous cosmologies, has culminated in the BKL- 
scenario Q. According to this scenario, as the singular- 
ity is approached, the spatial geometry can be viewed 
as a collection of small patches each of which evolves es- 
sentially independently as a homogeneous model, most 
generally the Bianchi IX model. This is justified by 
the observation that interactions between the patches are 
negligible because time derivatives dominate over space 
derivatives close to a singularity. 

The approach to the singularity of a Bianchi IX model 
is described by a particle moving in a potential with expo- 
nential walls (corresponding to the increasing curvature) 
bounding a triangle (Fig. During its evolution the 
particle is reflected at the walls resulting in an infinite 
number of oscillations (of the scale factors) when the sin- 
gularity is approached. This classical behavior can be 
shown to lead to a chaotic evolution by using an analogy 
with a billiard valid in the asymptotic limit close to the 
singularity |2|. 

To appreciate implications of the chaotic approach to 
the Bianchi IX singularity, observe that at any given time 
the spatial slice can be decomposed into a collection of al- 
most homogeneous patches, the size of the patches being 
controlled by the magnitude of space derivatives in the 
equations of motion. During subsequent evolution when 
the curvatures grow, these patches have to be subdivided 
to maintain the homogeneous approximation. This sub- 
division is also controlled by evolution of the individ- 
ual patches. Since the patches are homogeneous only to 
a certain approximation, a subdivision of a given patch 
at a certain time leads to new patches with slightly dif- 



ferent initial conditions. The chaotic approach to the 
Bianchi IX singularity then implies that their geometries 
will depart rapidly from each other, and the patches have 
to be fragmented more and more the closer one comes 
to the singularity. This rapid fragmentation suggests a 
very complicated and presumably fractal structure of the 
spatial geometry at the classical singularity. Note that 
both of these features, the breaking up of a spatial slice 
into approximately homogeneous patches and the unend- 
ing oscillatory approach to the singularity of individual 
patches, are ultimately consequences of the unbounded 
growth of the spatial curvature, i.e. the singularity. 

However, the classical evolution towards the singular- 
ity is expected to break down when the curvatures be- 




FIG. 1: The potential is shown for V = 30 in Planck units. 
Along the z-axis is plotted the logarithm of the potential 
(shifted so that it is larger than 1 everywhere). Its form 
does not change with volume due to the factorized volume 
dependence. 
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come too large and to be replaced by quantum dynam- 
ics. If we truncate the classical model at a certain lowest 
volume, technically chaos does not occur. It is a nat- 
ural question to ask if the qualitative features of the 
chaotic approach continue to survive quantum modifi- 
cations. The answer to this question is neither obvious 
nor independent of concrete theories of quantum grav- 
ity and must be explored within the context of a specific 
candidate theory. 

For instance, effective actions with extra terms and 
additional fields motivated from string or M theories do 
not lead to a non-chaotic behavior [5j and thus retain 
their implications for fragmentation in the context of a 
general inhomogeneous singularity. On the other hand, 
loop quantum gravity 0,0] predicts a discrete structure 
of space. For these theories, an unlimited fragmentation 
of space would be inconsistent with the discrete struc- 
ture and one would expect a qualitative modification of 
the chaotic approach. This translates into a self consis- 
tency test for such theories. Thus, a quantum theory of 
gravity with a discrete structure of space must, for self- 
consistency, provide a mechanism which prevents the un- 
limited fragmentation in the expected general approach to 
a classical singularity. 

In this paper we use methods of loop quantum cosmol- 
ogy 000- a part of loop quantum gravity, to study this 
issue. This allows us to obtain explicit, non-perturbative 
modifications of the classical behavior at small volume 
which can be analyzed for their implications for chaos. 
Loop quantum cosmology has already lead to a resolu- 
tion of conceptual problems such as a non-singular evo- 
lution and also given a new scenario for inflation [Io| 
which is based on the small-volume modifications. The 
proof of absence of singularities extends to homogeneous 
models, in particular the Bianchi IX model 00. The in- 
vestigations of this paper will provide a new consistency 
check of loop effects and their physical viability. Specifi- 
cally, quantum modifications for the vacuum Bianchi IX 
model (central to the issue of a chaotic approach) are 
presented and shown to prevent the chaotic behavior. 

The dynamics of the Bianchi IX model can be formu- 
lated on its minisuperspace spanned by the positive scale 
factors a/, / — 1, . . . , 3, related to the diagonal metric 
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It is given by the Hamiltonian 



H = - 
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where k = 8irG is the gravitational constant and 

r 1 f aj a K a] \ 

I> = 1 — for eij K = 1 (2) 

2 \a K a j cljclk / 

are the spin connection components. Thus, the potential 
term obtained from is given by, 



W{a l7 a 2 ,a 3 ) = ~ 



In order to diagonalize the kinetic term of QJ one 
can introduce the Misner variables the loga- 

rithmic volume £1 := — ^\ogV = — | \og(aia203) 
and the anisotropies 0±, introduced via a± := 

e -n+p + +^P-^ a2 .- e -n+p + -V3P^ a3 := e -n-2/3 + _ In 

terms of these variables the potential © is obtained as, 



W(Q,/3 + ,/?_ 



1 



4e- 2/3 +cosh(2V3/?_) + 
2e 4/3 + |cosh(4\/3/3_) - 1 j 



(4) 



(3) 



Evidently, the ^-dependence factorizes and one obtains 
an anisotropy potential V (/?+,/?-) which exhibits expo- 
nential walls for large anisotropies as displayed in Fig. ^ 
A typical wall can be derived from the potential by 
setting /3_ — and taking (3 + to be negative, e.g. 
W{a x ,ai,V/a\) = I e -^( e -8/3+_4e^ 2 ^) ~ ie" 40 " 8 ^. 

The wall picture implies that the universe, described 
by a particle moving in such a potential, runs through 
almost free (Kasner) epochs where the potential can be 
ignored, interrupted by reflections at the walls where 
the expansion/contraction behavior of different direc- 
tions changes. The infinite number of these reflections 
implies that the system behaves chaotically. 

One can see that the infinite height of the walls is a 
consequence of the diverging intrinsic curvature in Tj for 
small and large a/. In the classical evolution each a/ will 
eventually become arbitrarily small as the singularity is 
approached, e.g. a\ and a 2 for the typical wall above. 
However, at a certain stage of the evolution quantum 
gravity is expected to become important which will lead 
to a modification of the behavior. Clearly, a necessary 
condition for the chaotic behaviour to be prevented by 
the quantum modification is that quantum gravity should 
effectively contain an upper limit on the curvature such 
that the walls would have only finite height, changing the 
whole scenario. 

A quantum theory where such a maximal curvature 
follows is loop quantum cosmology [l2l Il3| . The origin 
of this upper bound for the curvature lies in the quan- 
tization of the relevant quantities and can be compared 
conceptually to the finite ground state energy of the hy- 
drogen atom obtained after quantization. In the classical 
equations of motion one can incorporate this feature of 
loop quantum cosmology by replacing the T/ with effec- 
tive coefficients which are derived from the quantization 
0. This leads to the corresponding effective potential. 

The new, effective potential is more complicated than 
the original one. Nonetheless, for large volume, the effec- 
tive potential approximates the original one. It is signif- 
icantly different for small volumes and is responsible for 
the breakdown of chaos. The volume dependence does 
not factorize, making the analysis of the classical mo- 
tion (with volume as internal time) harder. One should 
also keep in mind that the analysis done here uses only 
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the effective classical description which includes some 
non-perturbative quantum effects in the potential. The 
full quantum evolution is much more complicated and is 
given by a partial difference equation for the wave func- 
tion 0. 

Using the effective potential in classical equations of 
motion is sufficient to shed light on the above consistency 
requirement, even though it is not valid close to the clas- 
sical singularity. The underlying dynamical equation of 
loop quantum cosmology contains a parameter, j, that 
appears as a quantization ambiguity 14;]. Its value con- 
trols the size of the universe where the maximal curvature 
is attained. By choosing it to be sufficiently large, one 
can move the quantum effects in the effective potential 
into the semiclassical domain. Those large values may 
not be expected or natural from a physical point of view, 
but they allow us to study the consistency issue in a sim- 
plified setting. A necessary requirement for consistency 
of the quantum theory then is that the effective classi- 
cal description stops the Bianchi IX oscillations at some 
volume. 

To be more explicit, one first has to introduce den- 
sitized triad variables = V/aj, with the volume 
V = 0,1(120,3 = y / \p 1 P 2 P 3 \, which become basic opera- 
tors in loop quantum gravity. For convenience, we have 
taken them to be dimensionless by setting = 1 in 

the notation of 8]. Their inverses Ip 7 ! -1 do have well- 
defined quantizations as finite operators, despite the clas- 
sical curvature divergence at p 1 = which is eliminated 
by quantum effects 12]. One can model its effect in the 
classical behavior by replacing (p 7 ) -1 in the spin con- 
nection components by a function F(p I /2j) to get an ef- 
fective spin connection which includes non-perturbative 
quantum effects. The function F is derived directly from 
a quantization [Tij ] . The ambiguity parameter j appears 
explicitly and controls the peak of F. 

The potential W then gets replaced by an effective one, 
Wj, by using the effective instead of the classical spin 
connection components. One can again obtain the typical 
wall behavior, this time by evaluating Wj(p 1 ,p 1 , (V/p 1 ) 2 ) 
as a function of p 1 at fixed volume. At these arguments, 
with p 1 large enough but V/p 1 not too small, the compli- 
cated expression simplifies to just two dominant terms. 
As a function of p 1 and the volume, the typical wall be- 
comes (see Fig. 

W 3 {p\p\2 ]q ) « ^J^{3- 2qF(q)}, (5) 

where we have used the notation q :— ^{V/p 1 ) 2 . 

For values at the peak or larger, F(q) ~ g -1 , and 
the classical wall ie~ 4 °~ 8/3 + is reproduced. The peak 
of the finite walls is reached for a constant argument q 
of F which in Misner variables implies that e - 2 °+ 2 ^+ 
is constant. Thus, the wall maxima lie on the line 
f3 + = il + const in the classical phase space, and the 
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FIG. 2: Effective wall © of finite height for j = 1, V = 10 
as a function of x = — f3+, compared to the classical exponen- 
tial wall (upper dashed curve). Also shown is the exact wall 
Wj(p 1 ,p 1 , (V/p 1 ) 2 ) (lower dashed curve), which for x smaller 
than the peak value coincides well with the approximation © 
up to a small, nearly constant shift. 

height of the wall drops off as e _12n cx V 4 with decreas- 
ing volume. At very small volume, however, the walls col- 
lapse even more rapidly, and a numerical analysis shows 
that the potential becomes negative everywhere at a di- 
mensionless volume of about (2.172j) 3 / 2 in Planck units, 
i.e. just around the elementary discrete volume for the 
smallest value j = 1/2. 

In Fig. |21 are shown three snapshots of the effective 
potential at decreasing volumes in the vicinity of the 
isotropy point in the anisotropy plane. The potential at 
larger volumes clearly exhibits a wall (positive potential) 
of finite height and finite extent. As the volume is de- 
creased the wall moves inwards and its height decreases. 
Progressively the wall disappears completely making the 
potential negative everywhere. Eventually, the potential 
approaches zero from below. 

This immediately shows that the classical reflections 
will stop after a finite amount of time, rendering the 
classical argument about chaos inapplicable. The vol- 
ume where the transition from the classical behavior to 
the modified behavior takes place, i.e. the first time the 
universe can "jump over the wall," depends on the initial 
conditions but it will certainly happen — the latest at 
the elementary discrete volume for the smallest value of 
the ambiguity parameter. 

Starting at large volumes, the universe will pass 
through several phases of Kasner evolution punctuated 
by reflections. This will continue until the small volume 
modification of the effective potential comes into play. 
The stability of a Kasner trajectory in the presence of the 
modified potential will decide the nature of subsequent 
deviations after the 'last' Kasner epoch. Pleasingly, the 
Kasner trajectories turn out to be stable in the presence 
of the effective potential By contrast, these are un- 
stable in the presence of the classical potential. 
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FIG. 3: The modified effective potential is plotted for j = 5 
at volumes 150, 100, 30 in Planck units for —1.5 < {3± < 1.5. 

Intuitively, this gives us a natural and consistent pic- 
ture of the approach to a classical singularity in a quan- 
tum theory with a discrete structure of space: At larger 
volume the system follows the complicated classical evo- 
lution with different Kasner phases interrupted by reflec- 
tions at the walls. As described before, in the context 
of an inhomogeneous space this leads to a fragmentation 



into smaller and smaller patches. Once quantum effects 
are taken into account, the reflections stop just when the 
volume of a given patch is about the size of a Planck 
volume, i.e. at the scale of the discreteness. Below that 
scale, a further fragmentation does not take place and 
the discrete structure is preserved. We again emphasize 
that the modifications used here become important at 
the latest when the Planck scale is reached implying con- 
sistency with the expectations from a discrete structure. 
Below this scale, though, the effective classical descrip- 
tion should be superseded by the quantum description. 
We also emphasize that the existence of a maximal cur- 
vature is a consequence of quantization and not put in, 
in an ad-hoc manner. The results thus constitute a con- 
sistency test for loop quantum gravity. 

This also indicates that the results of USUI) El which 
prove a non-singular quantum evolution of homogeneous 
models can be generalized to the full theory, removing 
also inhomogeneous classical singularities. At the present 
stage, however, the results for the general case are to 
be regarded as preliminary and have to be supported 
by more general techniques directly in inhomogeneous 
quantum models. 
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